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ABSTRACT 

' Type la Supernova Hubble residuals have been shown to correlate with host galaxy mass, imposing 

a major obstacle for their use in measuring dark ener gy properties. Here, we calibrate the fundamental 
O ■ metallicity relation (FMR) of lMannucci et alj (|2010D for host mass and star formation rates measured 

' from bro ad-band color s alon e. We apply the FMR to the large number of hosts from the SDSS-II 

sample of lGupta et al.l (j2011l ) and find that the scatter in the Hubble residuals is significantly reduced 
^ . when compared with using only stellar mass (or the mass-metallicity relation) as a fit parameter. 

Our calibration of the FMR is restricted to only star-forming galaxies and in the Hubble residual 
calculation we include only hosts with log(SFR) > —2. Our results strongly suggest that metallicity 
On , is the underlying source of the correlation between Hubble residuals and host galaxy mass. Since the 

FMR is nearly constant between z = 2 and the present, use of the FMR along with light curve width 
and color should provide a robust distance measurement method that minimizes systematic errors. 
Subject headings: supernovae: general 



O 

U 



o 



00 



>< 



1. INTRODUCTION 



Oh. Type la supernovae (SNe la) are important for their use as distance indicators. They are not perfect s tandard 
cand les, however; the width of the light curve correlates wi th the peak brigh tness (i.e. the Phillips relation, iPhillipsI 
19931) . and color can be used to correct for dust extinction (jRiess et al.lll996l) . These correlations are instrumental in 
^ , the use of SNe la to measure the cosmolo gical properties o f our u niverse, and these mea surements have shown that 
our universe will not only expand fore ver (iGarnavich e t alT 1998a; Pcrlmuttcr et al. 1998 ?), b ut th a.t the expans ion is 
— I actually accelera ting (jRiess et al.lll998l: iPerlmutter eTal.. 1999 : Knop et al.. 2003; Tonrv et 31.1120031 iBarris et al.]|2004l: 
iRiess et ani2004D . This acceleration, driven by a "dark energy" (so called because the source is unknown), has been 
one of the most surprising scientific discoveries of recent history. 



' Cosmological s tudies have shifted fro m confirming acc eleration to constrainin g the properties of dark energy 
; IGarnavich et al."1998b; Astier et al."2006; Ricss et al. 2007: Wo od- Vasev et al.l2007tlMikn aitis ct a l.l2007tlEisenstein et al.l 
H! ■ [2007: CQnley et al. 2011; S ullivan et a l. 2011: ,Suzuki_ct_al. 20j^ Yet as the number of SNe la has increased, the 



statistical error in the cosmological parameters has been decreas ed such that systematic errors have begun to dominate 
(|Wood- Vasev et al.l [20071: iKessler et aT]|2009al: iGuv et al.|[2010l ). The focus of SN la science has thereby shifted to 
CN . constraining other parameters that may affect their peak bright ness. One such pa r ameter that appears to be correlated 
' with Hubble residuals is the stellar mass of th e host galaxy (iKellv et al.l 120101 : iLampeitl et al.l 120101 : iSullivan et al.l 
[ [2010; Gupta et al. 2011). Sullivan et alJ' (|2006[ ) also showed that the SN la rate was dependent on the host galaxy 
stellar mass and star formation rate (SFR), indicating an environmental effect. While the underlying source of this 
correlation is unknow n, logic suggests tha t the stellar mass must be tracing some under lying property, such as the 
progenitor metallicity (|Timmes et al.|[2003l : iKasen et al.ll2009l) or age (jKrueger et al.|[20Tol) . 
Theory does predict th at the amount of Ni-56 produced may depend on the metallicity of the progenitor (|Timmes et al.l 
5^ ■ [200l IKasen et al.l [200l . in the sense that more metal-rich environments produce a higher fraction of stable Ni-58 
(due to excess available neutrons) and therefore less luminous SNe la. Since the light curve is driven by the decay of 
Ni-56, an increase in the amount of stable Ni-58 produced will decrease the amount of Ni-56, and should reduce the 
peak brightness of the SN and change the shape of the light curve. Most of this correlation should be accounted for by 
the Phillips r elation, but with s o me residu al error due to small predicted variations in the Phillips relation at different 
mctallicities (jKasen et al.ll2009l ). IKasen et al. (20 0^ predict s that the systematic error on distance measurements due 
to this effect should be around 2%. iD'Andrea et al.l (|2011t l measured gas-phase metallicities of SN la host galaxies 
directly and found a strong c orrelation (Aa) w i th Hu bble residual^, suggesting that the stellar mass correlation could 
be a tracer of the metallicity. iGallagher et al.l (|2008l ) studied absorption line strengths in early-type host galaxies and 
found both an age and a metallicity corr elation with Hubble residuals. 

Recent work by iMannucci et ahl (|2010[ ). using Sloan Digital Sky Survey (SDSS) DR7 catalogs containing emission 
line fluxes and other spectroscopic properties of galaxief[3, discovered that galaxy metallicity is correlated not only to 
stellar mass, but also to SFR. Using the SFR in addition to the mass as a basis to estimate a galaxy's metallicity 
can significantly increase the precision of the metallicity estimate compared to mass alone (i.e. the Tremonti relation, 
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iTremonti efall 120041) . [M annucci et al.l (|2010[ ) found that star-forming galaxies up to z=2.5 foUow the same focal 
relation between metallicity, m ass, and SFR, wi t h no e vidence for evolution. Such a correlation, termed the fundamental 
metallicity relation (FMR) bv lMannucci et all ()2010[ ). can improve metallicity estimates when direct measurement of 
metallicity from relevant emission lines is not possible. While the FMR was o riginally defined for stellar masses a bove 
lO^-S M(^, it was later exte nded towards lower masses (jMannucci et al.|[2011| ). Among others. lYates et all (|2012[ ) and 
iPerez-Montero et al.l (120121 ) obtained qualitatively very similar relations between metallicity, mass and SFR, although 
with some quantitative differences. 

For stellar masses above about lO^"'^ Mq, the metallicity is nearly constant. Recent work bv I Johansson et all ()2012f ) 
found very little Hubble residual correlation with stellar mass for host galaxies with mass > 10^°'^ M©. This lack of 
correlation at higher masses could indicate that the mass is tracing the metallicity, considering the lack of variation 
with Hubble residuals in this region and the constancy of the metallicity. A significant goal of this paper is to test the 
ability of stellar population synthesis models to reproduce the FMR, which would then allow for accurate metallicity 
estimates using multi-band photometry alone. 

In this paper, we use the MPA/JHU catalogs of galaxy properties to calibrate the fundamental metallicity relation 
for use with t wo multi-band photometric stellar po pulation synthesis fitters. We present our me t hod a nd results 
using Z-PEG (jLe Borgne fc Rocca-Volmerangell2002t ). as well as the fittin g method of Gupt a et al.l ()2011[ ) using the 
Flexible Stellar Population Syn thesis (FSPS| models (jConrov et al.l 120091 : IConrov fc Gunnll2010( ). Hereafter, we wiU 
refer to the fits using the Gupta et al.l (|2011| ) fitter as Gil. The use of the FMR based on photometry alone allows 
accurate metallicity estimates of faint galaxies, or galaxies that are otherwise difficult to measure spectroscopically 
(for instance, high redshift galaxies where the relevant emission lines are redshifted beyond the capabilities of ground- 
based spectroscopy). This also allows for a study of the SDSS-H SNe la using the FMR, as the large sample of severa l 
hundred SNe la contains only about 40 galaxies that have emission-line estimated metallicities ('D'And rea et al.ll2011| ) . 
Using the improved metallicity estimates based on these photometric measurements, we test the Hubble residuals of 
the SDSS-H SNe la for their correlation with host metallicity. 

In section [21 we discuss the selection of the SDSS-H host galaxies, and the selection of the SDSS-H star-forming 
galaxies from the MPA/JHU catalogs. In section |31 we discuss the calibration of the FMR using Z-PEG and Gil, and 
report our results of testing the Hubble residuals against the metallicity parameter from the FMR. In section |4] we 
discuss some issues with the photometric fitters used in the paper, as well as extending the FMR to lower metallicities 
and higher redshifts. Finally, we present our conclusions in section [5l 

2. DATA 
2.1. SDSS host galaxy sample 

For this paper we have used the 206 SDSS-H host ga laxies from iGupta et al.l ()201lD . While a detailed summary 
of the exact method is presented in IGupta et al.l ()2011| ) , we will briefly summarize the SDSS-H SN Survey and the 
selection of the host galaxies in this section. 

Th e SDSS-H SN Survey searched 300 deg^ of sky, known as stripe 82, during 3 month periods of 2005 , 2006, and 
2007 (iFrieman et al . 2008) with a cadence of « 2 days using the CCD camera on the SDSS 2.5 m telescope ()York et al.l 



l2000t iGunn et al.l 1998, .2006i'). Around 500 tvp e la supernovae were spectroscopically confirmed over the course of the 
survey (jSako et al ![ l200ilHoltzman et al.|[20Cll . The SNe la range in redshift from 0.01 to 0.42, with a median around 



0.2. 

We use the publicly available Supe rnova Analysis package SNANA (jKessler et al.l l2009bl ). using the light-curve 
fitting process SALT2 (|Guv et aL 2007) t o determine t he best-fit distance modulus for each supernova from the SDSS- 



II photometry (jHoltzman et al.ll2008; Sako et al.ll2011|) . Initial cuts on the sample of SNe were made as follows: 



1. at least one measurement with Trest < —2 days, where Trest is defined as the number of days from the date of 
maximum; 

2. at least one measurement at Trest > +10 days; 

3. at least one measurement with signal-to-noise ratio (S/N) > 5 for each of the g, r and i bands; 

4. best-fit probability greater than 0.001, calculated from the value of the SALT2 fit. 

These cuts reduce the sample to 319 SNe. 

Each SN was visually inspected for host galaxy identification. 14 SNe were removed from the sample because they 
do not have identifiable hosts, either because they fall outside of the SDSS footprint, were too faint, or had r-band 
magnitudes outside the range 15.5 < r < 23. After ligh t-curve fittin g and host inspection the sample contains 305 
objects. The final cut is the best-fit probability from the IGupta et al.l (|2011[ ) fits. Fits with probability < 0.001 were 
removed, which leaves 206 SNe in the SDSS-H sample. 

2.2. SDSS general galaxy sample 

For comparison w ith iMannucci et al.l (I2010D . we have use d the MPA/JHU c atalogs for the SDSS-DR7, availa ble at 
http://www.mpa-garching.mpg.de/SDSS/ and described in iKauffmann et "ahl (|2003b[ ). iBrinchmann et all ()2004[ ). and 
.Salim et al.. (,2007.) . The catalogs contain 927,552 galaxies with measured spectroscopic properties, as well as stellar 
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masses (measured photometrically), multi-band photomet ry, and other importan t parameters. We have applied similar 
cuts to the data in order to reproduce the sample used bv lMannucci et aiT(|2010f) in the calculation of the fundamental 
metallicity relation (FMR). Our cuts are as follows: 

1. Redshift between 0.07 and 0.30. As described in iMannucci et al.l ()2010[) . this ensures that the [OII]A3727 line 
falls in the useful spectral range, as this line is important in calculating the metallicity, and that the 3-arcsec 
aperture of the spectroscopic fiber samples a significant fraction of each galaxy. 

2. The signal-to- noise ratio of Ha is > 25. This requirement leads to a more reliable value of the metallicity. In 
the catalogs it is mentioned that all Ha errors can optionally be scaled by 2.473; we do not perform this error 
modification as it cuts many more galaxies. If this error scaling is correct, this S/N cut on Ha would be roughly 
10, instead of 25. 

3. The Balmer decrement must be greater than 2.5, and the dust extinction. Ay, must be less than 2.5. The 
extinction was calculated from the Balmer decrement as: 

Av = 3.1 X 2.1575 x logio ( ^""^Se""' ) 

4. AGN were removed using the criteria in IKauffmann et "all (|2003aD . 

5. Beyond the IMannucci et all ()2010l ) selection criteria, we also required that 

[OII]A3726 + [OII]A3729 + [OIII]A4959 + [OIII]A5007 > 0, 

and Ha and H/3 > 0. The sum of these oxygen lines is used in the calculation of the metallicity indicator R23. 
This ensures that no errors occur when we calculate the logarithm of these values in estimating the metallicity, 
as we had a small percentage of galaxies which did not pass these cuts. 

We calculate the SFR from the Ha emission-line flux, corrected for dust extinction estimated from the Balmer decre- 
ment. We use the iKenn icutt' (1998) conversion factor between Ha luminosity and SFR. While IMannucci et al.l ()2010f ) 
cor rects this value from a Kroupa (2001) IMF to a Chabricr (2003) IMF, wc do not modify the value because we use 
the lKroupal (pOfH IMF in our SED fitting process, following previous SN la host galaxy analyses. 

For the oxygen gas-phase abundanc e we use a comb i natio n of R23 = ([OII]A3727 -I- [OIII]A4958,5007)/H/? and the 
[NII]A6584/Ha ratio, as described in iMaiolino et al.l (|2008l ). We restrict the values such that log([NII] A6584/Ha) 
< —0.35 and log(R23) < 0.9, and then select only galaxies where these two abundances are within 0.25 dex of each 
other. These restrictions resolve the ambiguity in the R23 metallicity estimate. 

Our final sample contains 142,227 galaxies, as compared to 141,825 galaxies in the IMannucci et aLl ()2010l ) sample. 

2.3. Photometric fitting using Z-PEG and FSPS Templates 

For photometric fitting of ugriz photometry of galaxies, we use the code Z- PEGF| ()Le Borgne fc Rocca-Volmerangd 
I2002D . as weU as the FSPS model minimization technique of iGupta et"all (|201ll ) (Gil). We use version 5.23 of Z- 
PEG, obtained via private communication from Damien Le Borgn^. Z-PEG utilizes the PEGASE.2 stellar population 
synthesis code, using 9 different star formation scenarios to build models of galaxy spectral types from elliptical to 
starburst. Synthetic spectra are produced at galaxy ages of to 20 Gyr. The code performs a minimization over 
all templates convolved with the SDSS filter functions in order to select the best-fit galaxy model, and the parameters 
of this model such as stellar mass and SFR are considered the be st-fit parameters of the galaxy. Error s on these 
parameters are calculated from the information of the full space ()Le Borgne fc Ro cca-Volmcrangc"'2 002D. 

For our fits, we generally use the default parameters as described in table 1 of [Le Borgne & Roc ca-Volmerang^ 
(gOO^. There are two main differences; we leave out the starburst template, and we average the SFR over the past 
500 million years. Both of these differe nces are com monly ap plied in the litera.ture f or using Z-PEG to fit type la 
supernova host galaxies (iSuUivan et al.l [2 006: Lanipe itl et al.l [2 010: S mith et al.ll20lj] ). We use the other 8 galaxy 
evolution scenarios described in lLe Borgne fc Rocca-Volmeranga ((2002i) . with 200 age steps instead of 69, in order to 
provide better age resolution for older galaxies (the 69 time step scenarios have a time difference of 1 Gyr after ages of 
about 5 Gyr). For stellar mass fitting, we utilize both redshift constraints; the first constraint is that the galaxy cannot 
be older than the age of the universe at the given redshift, and the second is that the galaxy must be at least 11 Gyr 
old at z = 0. The second constraint is a little more relaxed than that listed in lLe Borgne fc Rocca-Volmerangd ((2002), 
but still has the effect of requiring the best-fit template to contain old stars. For fitting of the SFR, we remove the 
age constraint at z=0. This allows templates with higher specific star-formation rates, and increases the diversity of 
the best-fit templates in sSFR. Testing indicated that this method best matched the SFR and sSFR as measured from 
emission lines. We also use the same cosmology as Gupta et al. (20 0), with Q\ = 0.735, — 0.274, and Hp = 0.70. 
This c osmology is taken from the SDSS-II SN Survey year 1 analvsis ()Kessler et al.ll2009aD . Lastly, we use the iKroupal 
(pool initial mass function, following the method of the previously mentioned type la supernov a host analyse s. 

For the Gil fits, the models were built using the same FSPS parameters as listed in table 1 of IGupta et al.l (|2011| ). 
A x'^ minimization is performed to select the best-fit model. 

http://imacdlb.iap.fr:8080/cgi-bin/zpcg/zpeg.pl 

® Information about PEGASE.2, Z-PEG, and contacting tiie developers can be found at |http: //www2.iap.fr/pegase/| 
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Figure 1. The fundamental metallicity relation using the stellar masses from the MPA/JHU catalogs. Stellar masses were photometrically 
estimated, and the SFR is calculated from the Ha emission. Left: the metallicity of the SDSS-DR7 galaxies, as measured from the 
MPA/JHU catalog emission lines, as a function of the catalog photometric estimate of the galaxy mass. The points represent the median 
of the metallicity and the mean of the mass, in bins 0.15 dex wide in the mass. Right: The same objects, also with 0.15 dex wide bins in 
the emission-line-calculated SFR, showing the reduction in scatter that the SFR provides. 



3. ANALYSIS 

Our ultimate goal is to study the correlation of SN la Hubble residuals with host galaxy stellar masses, using the 
fundamental metallicity relation, in order to determine the underlying source of this correlation. We will calibrate 
the FMR using masses and SFRs from two multi-band SED fitting procedures; we use multiple fitters to show the 
robustness of the FMR using different stellar population models. The use of the FMR in our study is necessary because 
our sample of host galaxies contains a small fraction of spectroscopically measured metallicities. We will calculate the 
FMR for each of these fitters, and then independently estimate the effect that both the metallicity converted from the 
FMR, and the FMR parameter, /xfmr — log(mass) — ax log(SFR), have on the Hubble residuals of our sample. 



3.1. Reproducing the Mannucci FMR 

We begin by reproducing the FMR of iMannucci et aP (|201Cll ). The primary reason for reproducing this analysis is 
to set the framework and the basic comparison for all photometric formulations of the FMR. For this analysis we use 
the stellar masses from the MPA/JHU catalogs, and we use the SFR and metallicities as measured from the emission 
lines in the catalogs, described in section [2^ We bin the 142,227 galaxies in 0.15 dex wide bins in both log(mass) 
and log(SFR). The metallicity in the bin is estimated by calculating the median value. For the cent er of each bin in 
stellar mass and SFR, we use the mean of the values within the bin. This is slightly different from IMannucci et al.l 
(|2010( ) because we found that the center of the bin was not an accurate representation of the center of the data for bins 
that were near the edge of the distribution. For the error in metallicity in each bin, we use the standard error of the 
median, given by Se = : where cr is the standard deviation of the metallicity values in the bin and N is the number 
of galaxies in the bin. slight differences in building the sample and in the binning process require a re-calibration of 
the FMR in order to set the basis for comparison with photometric techniques. 

We then fit a 4th order polynomial to the median meta llicity versus mass dis tribution. This is shown in the left pane 
of figured] This is analogous to the Tremonti relation (jTremonti et al.l [20041 ). This 4th order polynomial fit has an 
RMS dispersion of 0.065 dex. We next fit a 4th order polynomial to the median met allicity versus Mfmr dis tribution, 
where /ifmr = log(mass) — ax log(SFR). We use the label /ipMR instead of /i from IMannucci et al.l ([2010), because 
in SN la cosmology the parameter /i usually represents the distance modulus. For all FMR and mass-metallicity 
relation fits reported in this paper, we minimize the value of a 4th order polynomial using the Nelder-Mead simplex 
algorithm minimization in the Pythor|3 module scipv. optimiz J^. allowing the FMR parameter a to vary in the fit. 
The best fit 4th order polynomial and a value are given by: 

12 -I- Log(0/H) = 8.92 + 0.365x - 0.246x2 - O.OSlx^ -I- O.lOQx'' (2) 

with X = hfmk — 10, and ^fmr, = log(mass) — 0.346 x log(SFR). This best-fit polynomial is shown in the right pane 
of figure [U and reduces the RMS dispersion from 0.065 to 0.02 dex. This reduction also appears significant because 
the metallicity estimate appears to be more gau s sian a nd symmetric about the best-fit line. Our measurement of the 
FMR matches well with that of Mannuc ci et al.l (|2010[). wh o found a — 0.32 with a minimum dispersion of 0.02 dex. 

We have reproduced the results of Mannucci et al.l (|2010D . in showing that adding a fraction of the SFR can reduce 
the scatter in metallicity estimates of galaxies. This formulation of the FMR sets the basis for comparisons with the 
photometrically estimated FMR. 

® http://www.python.org/ 

http:/ /docs. scipy.org/doc/scipy/reference/generated/scipy.optimize. fmin.html#scipy.optimize.fmin 



5 




3.78 




3.36 




2.94„ 

tn 




V 

2.52 a 

§1 


zpeg 


2. ID'S 

V 

1.68 1 

3 
C 

1.26 0; 




log(SFR) 


0.84 




0.42 




0.00 





8 9 10 11 

log(mass) SDSS catalog 




-2-10 1 2 
log(SFR) SDSS catalog 



Figure 2. Density plots (greyscale levels in log space) of the Z-PEG fits compared to the catalog mass and the SFR measured from 
Hq. Left: the Z-PEG stellar mass estimates against the SDSS DR-7 mass estimates from'Sal im et all 112007 1). The offset is due to the age 
constraint at z = used in Z-PEG. Right: the Z-PEG SFR values against the values measured directly from emission lines. Z-PEG 
templates with no recent star-formation were assigned an arbitrary value of log(SFR) = —2.5 in this figure. These objects are not included 
in the calculation of the FMR. 
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Figure 3. Density plots (greyscale levels in log space) of the Gil fits compare d to the catalog m ass and the SFR measured from Hq. 
Left: the Gil stellar mass estimates against the SDSS DR-7 mass estimates from lSalim et all 120071) . There is also a slight offset in these 
values, possibly due to the use of a different IMF. Right: the Gil SFR values against the values measured directly from Hq emission lines. 



3.2. Reproducing the FMR with Photometric Fitting Techniques 

In figure [21 we compare the mass and SFR derived from the Z-PEG fits to the SDSS DR-7 sample. The mass values 
are highly correlated but have a sy stematic ofi^set of .20 dex as compared to the catalog, which were also calculated 
photometrically by the method of iSalim et al.l (|2007| ). This offset is caused by the age constraint at z = in the 
Z-PEG fits. Relaxing the age constraint brings the main group of galaxies closer to the catalog masses, but has the 
effect of increasing the number of outliers at preferentially low masses. The SFRs as calculated by Z-PEG also show 
a systematic offset with respect to the SFR measured from Ha. Z-PEG templates with no recent star formation 
were assigned an arbitrary value of —2.5 for the contour plot. These objects are removed from the sample prior to 
calculation of the FMR, along with objects with poorly constrained stellar masses, reducing the sample size to 140,987. 

Figure |3] shows the mass and SFR values from the Gil fitting method against those from the catalog. While the 
mass values show a small offset similar to the offset from Z-PEG, this offset is likely due to the use of a different IMF 
in the Gil fitting method. The SFR values show a tendency to have less spread, and with a less-steep slope. Slightly 
different from Z-PEG, Gil does not generate zero SFR objects. The objects that extend down to low SFR values are 
likely the same objects that result in zero SFR galaxies in the Z-PEG fits. There is also a curious group of galaxies 
in stellar mass space, with much larger values than the catalog. However, the number of galaxies in this region is 
relatively small compared to the full sample, so this does not systematically affect the FMR. The greyscale levels are 
drawn in log space, so the number of galaxies in this region is very small compared to the main sample. 

We use the same binning method as in section 13.11 to calculate the median metallicity values as measured from the 
emission lines. Using Z-PEG to fit the ugriz multi-band photometry for the 140,987 SDSS galaxies passing our cuts, 
we reproduce the reduction in scatter as a result of combining mass and SFR to estimate metallicity. Figure |4] shows 
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Figure 4. The fundamental metallicity relation using Z-PEG applied to broadband photometry. Adding a portion of the SFR significantly 
reduces the scatter in the metallicity estimate. All mass and SFR values were estimated photometrically by Z-PEG, from ugriz photometry 
from SDSS-DR7. This demonstrates the power of photometric fitting of multi-band photometry to estimate the metallicity of a large 
sample of galaxies. Left: the metallicity of the SDSS-DR7 galaxies as a function of the Z-PEG estimate of the galaxy stellar mass. The 
points represent the median of the metallicity and the mean of the mass, in bins 0.15 dex wide in the mass. Right: The same objects, also 
with 0.15 dex wide bins in the Z-PEG estimated SFR, showing the reduction in scatter that the SFR provides. 




Figure 5. The fundamental metallicity relation using the Gil fitting method. Adding a portion of the SFR, significantly reduces the scatter 
in the metallicity estimate. All mass and SFR values were estimated photometrically using FSPS stellar models in the minimization 
routine of Gil, from ugriz photometry from SDSS-DR7. This demonstrates the power of photometric fitting of multi-band photometry to 
estimate the metallicity of a large sample of galaxies, and corroborates our findings using Z-PEG. Left: the metallicity of the SDSS-DR7 
galaxies as a function of the Gil estimate of the galaxy mass. The points represent the median of the metallicity and the mean of the 
mass, in bins 0.15 dex wide in the mass. Right: The same objects, also with 0.15 dex wide bins in the Gil estimated SFR, showing the 
reduction in scatter that the SFR provides. 



the median metallicity as a function of mass (i.e. the Tremonti relation of iTremonti et all 120041 ) . and as a function 
of log(mass) — 0.306 x log(SFR). The FMR reduces the dispersion from 0.05 in the mass-only case to 0.019. The 
equation of the best-fit line to the median metallicities is: 

12 + Log(0/H) 8.94 + 0.284x - 0.226x2 - O.OlOx^ + 0.058x* (3) 

with X = ^FMR - 10, and fipMR — log(mass) — 0.306 x log(SFR). 

Figure [S] shows the correlation of Log(mass), and /ipMR = Lo g(mass) — 0.534 x Lok( SFR) for Gil fits. Once again, 
the dispersion is significantly reduced, confirming the results of IMannucci eFall (pOlOf ) and agreeing with the Z-PEG 
results in the sense that photometric fitting of mass and SFR values can reproduce the FMR. The stellar mass-only 
metallicity estimates have an RMS value of 0.065 dex, while the FMR estimated metallicity values have an RMS 
dispersion of 0.03 dex. The Gil FMR is given by: 

12 -I- Log(0/H) = 8.81 -I- 0.718X - 0.598x2 _ o.074x3 -|- 0.223x'* (4) 

with X = fiFMR - 10, and /xfmr — log(mass) - 0.534 x log(SFR). 

Figure m shows the dispersion in the metallicity as a function of the factor a, where the fit parameter for the FMR 
is given by log(mass) — axlog(SFR). The minimum value occurs where a = 0.306 for Z-PEG , and a = 0.534 for 
Gil fits. There is a difference in the FMR as calculated from Gil fits from IMannucci et al.l (|2010D as well as our own 
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Figure 6. The per degree of freedom value of the best-fit polynomial as a function of the parameter a, from the Z-PEG FMR (left), 
and the Gil FMR (right) calibrations. The metallicity estimate is a function of log(mass) - axlog(SFR). For this plot we fix a at the 
indicated x-axis value and calculate the best fit fourth-order polynomial. 

analysis of the emission lines, which results in a arou nd 0.32. We discuss this issue in section S) 

We reproduce the results of 'Man nucci et all (j201C(l . using photometric measurements alone, in that using both the 
st ellar mass and a fractio n of the SFR reduces the scatter in the metallicity estimates over the mass- metallicity relation 
of iTremonti et all ()2004[) . The mass-only metallicity estimates for Z-PEG have an RMS dispersion of 0.05 dex, while 
the FMR estimates have a dispersion of 0.019 dex. For FSPS templates in the fitting method of Gil, the FMR reduces 
the dispersion from 0. 06 to 0.03. This is compared to the RMS dispersion using the SDSS catalog of 0.02 dex, as in 
IMannucci et al.l (|2010( ). We have shown that photometrically estimated SFRs are sufficient to significantly reduce the 
scatter in the metallicity estimates using only multi-band photometry. 

3.3. Applying the FMR to Hubble Residuals 

In this section we discuss the application of the FMR to the Hubble residuals of the SDSS-II SNe la in our sample. 
We present the analysis for our Z-PEG fits, and in sectionUwe will present the same analysis for Gil fits. As a basic 
outline, we use two primary methods to compare the mass correlation with the FMR. First, we use both the mass- 
metallicity relation and the FMR to convert the SDSS-II host galaxy photometric measurements to metallicity. We then 
determine the correlation between these metallicities and the Hubble residuals, and compare them. Then, we compare 
the Hubble residual correlation of mass alone with the Hubble residual correlation of ^fmr = Log(mass) — ax Log(SFR) 
, allowing a to vary. 

To fit a line to t he Hubble re siduals as a function of any parameter, we use the Monte Garlo Markov Ghain (MGMG) 
method LINMIX (|Kellvl [20071 ). This allows for a robust measurement of the linear fit parameters while taking into 
account the error in the independent variable. Similar to iGupta et al.l ()2011f ). we use 100,000 MGMG realizations, 
and LINMIX provides the posterior distribution of the slope and intercept values of the best-fit line. Using these 
distributions, which are highly gaussian, we fit a normal distribution to each of the resulting posterior distributions 
for the slope and intercept. The mean of the gaussian is considered the best-fit value (for either the slope or the 
intercept), and the a of the gaussian is considered the error. Dividing the mean by this cr gives an indication of the 
statistical significance from zero slope, in terms of the number of standard deviations. 

While the mean and a values of our gaussian fits determine the statistical significance of any potential non-zero 
slope, this is not a quantitative measurement of the goodness-of-fit. To estimate how well the line fits the data we use 
the statistic. This allows us to directly compare multiple different independent variables based on the resulting 
value. In calculating the value of the best-fit line, we use both dependent and independent error sources; for the 
independent variable errors, we propagate the error by fixing the value of the FMR parameter a at the value of the 
best-fit FMR. As part of our analysis, we perform statistical testing as a function of a and fixing the value for error 
calculations prevents the error bars from varying with this parameter. 

When comparing a Hubble residual fit of the mass alone to a fit using a version of /^fmRj we use the likelihood ratio 
test (LRT). This statistic is used for nested models, in which the alternative model can be transformed into the base 
model by a simple restriction on the parameter space. The test statistic for the likelihood-ratio test. A, is given by: 



A = -21n 



(5) 



where Pq is the probability of the best-fit base model and Pa is the probability of the best-fit alternative model. This 
test statistic obeys the distribution, with degrees of freedom k = dfa — dfo, where df g and dfp are the degrees 
of freedom of the alternative and base models, respectively. The Neyman-Pearson lemma ([Nevman fc Pearsonlll933l) 
indicates that this statistic has the highest power when distinguishing between two models. When we report a likelihood 
ratio test probability, it is the result of using this parameter A as the test statistic in the x^ distribution. Usage of 
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Figure 7. The Hubble residuals as a function of metallicity converted from mass (left), and the fundamental metallicity relation (right) 
for Z-PEG. The blue line is the linear fit to the Hubble residuals, and the red dots are the median values of the Hubble residuals. The 
likelihood-ratio test shows that the value of the best fit line with the FMR metallicities is significantly better than that of the mass-only 
values, at a significance of 99%. 



this statistic in comparing the mass alone with the FMR is a test of the effect of the FMR parameter a on the Hubble 
residuals. In this scenario, a LRT probability less than 5% indicates a statistical preference for the alternative model 
(the FMR). In all of our tests we use a version of the FMR as the alternative model, and the mass alone as the base 
model. The base model is a simple restriction of the alternative model in that mass alone is simply the FMR parameter 
Mfmr with a = 0. If a and therefore the FMR does not improve the goodness-of-fit of the best-fit line to the Hubble 
residuals, the LRT will return a value that is not statistically significant. 

Using our Z-PEG fits, we remove host galaxies with log(SFR) < -1.8, leaving 193 galaxies from the original 206. 
Then we remove galaxies which have /ifmr values that fall outside the region of our photometrically calibrated FMR. 
This leaves 182 galaxies in the Z-PEG sample. For these galaxies left in our sample, we convert the photometrically 
measured stellar mass and SFR values using both the Z-PEG mass-metallicity relation and Z-PEG FMR. We then 
used LINMIX to fit the Hubble residuals as a function of both of these metallicities. In this analysis, we do not include 
the traditional 0.14 magnitude systematic error in the Hubble residuals. This error is intended to account for unknown 
sources of syste matic error, and it is an unknown source of systematic error that we are testing. This reasoning is 
similar to that of iGupta et al.l (|2011| ). 

Using the mass-metallicity relation values, we obtain a best-fit line with value of 481.08 with 180 degrees of 
freedom. This best- fit line has a non-zero slope with significance of 3.18(7. Using the FMR metallicity values, we 
obtain a best-fit line with value of 468.5 with 179 degrees of freedom. This non-zero slope has a significance of 
3.10(7. Using the likelihood ratio test to compare these fits, with the FMR parameter a as the lone degree of freedom, 
produces a LRT statistic A of 6.98 with one degree of freedom. In the distribution this results in a probability value 
of 0.008, indicating that the improvement in due to the addition of the FMR parameter a is significant at greater 
than 99% confidence. 

The fit of the Hubble residuals versus mass-metallicity relation has a slope of —0.222 (0.070) with an intercept of 
1.966(0.620). The fit with respect to the FMR metalhcity has a slope of -0.271 (0.087) with an intercept of 2.420 
(0.782). Figure [7] shows the Hubble residuals versus the metallicity estimates for both scenarios. The metallicity 
reaches a maximum near 9.07, and many galaxies are grouped in this region. The inclusion of part of the SFR in 
the metallicity calculation tends to move points left to right in these plots, and the LRT indicates that the fit to the 
FMR values is statistically preferred. At first glance, it seems odd that two fits that both have a 3a significance of 
non-zero slope would produce a LRT probability that prefers one over the other. It is important to remember that 
the significance of the non-zero slope is in no way an estimate of the goodness-of-fit. It only contains information 
about whether or not there is a correlation between the dependent and independent variables being tested. It cannot 
necessarily provide information about which independent variable is better at reducing the scatter in the distribution. 

The significance of our result is two-fold. First, we have shown that the Hubble residuals of 182 SDSS-II are strongly 
correlated with the metallicity estimates from both the mass and the FMR. We have also shown that the FMR provides 
a statistically better fit to the Hubble residuals, an indication that the stellar mass correlation with Hubble residuals 
is a tracer of a metallicity correlation. 

Because galaxies in the sample have a maximum metallicity value which causes galaxies to bunch up at high 
metallicity, we also test the Hubble residuals against the mass values, and against the /ipMR — Log(mass)— axLog(SFR) 
values over a range of a values. Figure [8] illustrates the LRT probability value as a function of a. To produce this 
plot, we fit the Hubble residuals versus hfmr for a values ranging from —0.1 to 1.5. We compared each fit to the 
case where a — using the value of the fit and the LRT. All errors were fixed to the a value for the best-fit FMR, 
simply to avoid the independent errors varying between fits. The fixed value of a in the error calculation does not 
affect these results. The vertical dotted line at a = represents the mass alone, and as values are compared to this 
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Figure 8. The likelihood ratio test probability as a function of the FMR parameter a, for Z-PEG. The probability value is a result of using 
the likelihood ratio test statistic A, from equation [5] as the test statistic in the probability distribution with 1 degree of freedom. This 
compares each value of a against the mass alone (shown by a vertical line at a=0). This minimization was performed using the Hubble 
residuals, and is mathematically independent of the FMR. The region of parameter space where the probability is below 5% agrees with 
the best-fit a value of the independently measured FMR (indicated by a star). This indicates that a linear combination of the stellar mass 
and SFR improves the Hubble residuals, and implies that the stellar mass correlation with SN la distances is tracing a correlation with the 
metallicity. 

value, the LRT probability for a = is 1. The vertical dotted line at a = 1 represents 1/sSFR, a parameter highly 
correlated with the age of the population. Probability values less than 5% indicate a significant improvement in the 
goodness-of-fit for the FMR parameter, suggesting a rejection of the null model (mass alone) in favor of the alternative 
model (non-zero a). The star indicates the a value of the best- fit FMR from the SDSS DR7 galaxies. The region 
of statistical preference in a agrees well with the independent formulation of the FMR, and does not favor the mass 
alone or the inverse of the specific star-formation rate. This is confirmation of the correlation of the Hubble residuals 
with metallicity, and further evidence that the Hubble residual correlation with stellar mass is tracing a correlation 
with the metallicity. 

We reiterate that since we are using the /ipMR value, log(mass) — ax log(SFR), the minimization that is performed 
to produce figure [8] via the likelihood ratio test is independent from the minimization performed when calculating the 
FMR. The only links are the log(mass) and log(SFR) values. The results show a definite improvement in the Hubble 
residual correlation when using Log(mass) and a fraction of Log(SFR). While not definitive, this strongly suggests 
that the Hubble residual correlation with stellar mass is actually tracing a correlation with the metallicity, in the case 
of star-forming host galaxies. 

4. DISCUSSION 

4.1. Comparing the Mannucci FMR to Z-PEG and \Guvta et al\ Ii201i ) 

We have demonstrated that the FMR can be formulated from broad band colors alone, but the parameters are 
dependent on the method used to estimate the stellar mass and SFR. In this section, we discuss some possible reasons 
for these differences in the functional form of the FMR. 

Figures IH [TOl and [11] contain some evidence as to the main difference between the estimated mass and SFR values 
that causes the fluctuation in a. These figures show the log(SFR) vs log(mass) for the catalog, Z-PEG, and Gil model 
minimization techniques. The simplest difference to notice on first glance is the systematically lower SFR values 
from Gil as compared to Z-PEG. However, testing indicated that applying systematic offsets to the SFR values did 
not significantly change the best-fit a parameter in the FMR. 

The next difference is the systematic offset in stellar mass values present in both Z-PEG and Gil as compared to 
the catalog. This effect can be tested by "mixing and matching" mass estimates in the FMR. In other words, using 
one measurement of the SFR with the three different measurements of the mass (catalog, Z-PEG, and Gil models). 
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Figure 9. Greyscale density plot (contour levels in log space) of the sample of galaxies with mass values from the SDSS DR-7 catalogs, 
and SFR values measured from Ha. 

Doing so results in a spread of a values around 0.1 for each SFR method, not nearly enough to explain the difference 
in the FMR compared to the catalog stellar masses and SFR from Ha. We conclude that the mass offset is not the 
main source of a variation between FMR calibration sets. 

The last major difference is in the specific star- formation rate. In figures [9l [TOl andfTTl we show the SFR versus the 
mass. The distribution is much narrower for the Gil fits than for the catalog or Z-PEG values, indicating a lack of 
spread in the specific star formation rate. The sSFR against mass for all datasets is shown in figure [121 The Gil fits 
show a marked difference in the sSFR distribution. Gil has no sSFR values above a certain level, which are clearly 
represented in the Ha measured values. This is caused by an age constraint in the Gil fitting process combined with 
a declining SFR model. Indeed the spread in sSFR values measured from Gil is around half that of the Ha measured 
values. This results in a significantly lower number of bins in the stellar mass / SFR plane, with the galaxies packed 
closer together in parameter space. This is the reason that the Gil FMR in figure[5]has many fewer points than figure 
[T] and figure ID In order to account for this mathematically, the minimization process compensates with a larger a 
value. Indeed "mixing and matching" the SFR with a fixed mass measurement method reproduces the large changes in 
a that we observe when using the different calibration methods to estimate the FMR. The small range in log(SFR) and 
log(sSFR) estimated from the photometric fits in Gil means that doubling a is required to get the same metallicity 
as measured directly from spectroscopy. 

4.2. \Guvta et al\ 1(201 i ) Likelihood Ratio Test Results 

In this section we briefiy reproduce the primary figures and results from the analysis section using Gil fits. Starting 
with the 206 host galaxies, removing galaxies with log(SFR) < -2 leaves 163 galaxies. Removing galaxies that do not 
fall in the calibrated FMR region further reduces the sample to 135 galaxies. The primary reason for the difference in 
the number of galaxies is the specific star formation rate distribution of the fits; this was discussed in section HTT] 

We first compare the best-fit value of the Hubble residuals versus mass-metallicity relation and the Hubble 
residuals versus FMR metallicity. The distributions can be seen in figure [131 The mass-metallicity best-fit line has a 
value of 356.84 with 133 degrees of freedom. The FMR metallicity best-fit line has a value of 353.08 with 132 
degrees of freedom. This results in a LRT statistic A = 1.39, with 1 degree of freedom. The probability of this value 
in the distribution is 0.239, a result that indicates no statistical preference for the FMR metallicity. 

The LRT as a function of the FMR parameter a using the mass and /j,fmr values once again agrees with the FMR, 
in the sense that including a portion of the SFR reduces the Hubble residual scatter at a statistically significant level. 
FigurellHindicates that adding an element of the SFR to the mass reduces the Hubble residual scatter about the best-fit 
line, for a values that agree well with the Gil formulation of the FMR. This again confirms the metallicity correlation 
with Hubble residuals, and suggests that the Hubble residual correlation with mass is tracing this correlation with 
metallicity. 

As discussed in the previous section, the Gil fits remove many more galaxies from the sample due to the lack 
of spread in sSFR values. This lack of spread in sSFR is what causes the LRT test for direct comparison between 
metallicity estimates to be statistically insignificant. The loss of around 50 galaxies from the sample simply due to 
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Figure 10. Greyscale density plot (contour levels in log space) of the sample of galaxies with mass and SFR values measured with Z-PEG. 
Shown in red are the SDSS-II hosts. The SFR distributions are different because of a bias towards high SFR galaxies due to our selection 
criteria in the SDSS DR7 catalog. 
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Figure 11. Greyscale density plot (contour levels in log space) of the sample of galaxies with mass and SFR values measured with FSPS. 
Shown in red are the SDSS-II hosts. 



their mass and SFR values falling outside of the FMR calibration range indicates the importance of improving the way 
SED fitters estimate the SFR and therefore the recent star formation history. Similar to the default Z-PEG setup, Gil 
uses an age constraint which forces each template to contain old stars. This forces an unrealistic correlation between 
the mass and SFR values. While this constraint reduces the number of outliers in stellar mass, but greatly restricts 
parameter space with respect to the recent SFH. 

For our Z-PEG fits, relaxing the age constraint was a simple process. Gi l, however, uses FSPS templates that do 
not have a published set of parameters that mimic real galaxies (similar to iLe Borgne fc Rocca-Volmerangell2002l for 
Z-PEG), so we were not able to reproduce the sSFR distribution of the catalog and Z-PEG values. Looking ahead 



12 




Figure 12. The specific star formation rate versus mass, for; left: the DR7 catalogs, middle: Z-PEG, and right: the fits to FSPS 
templates. The cutoff and lack of spread in sSFR values in the photometric fitting techniques explain the difference in the FMR a value 
between the catalog and Z-PEG FMR formulations, which agree quite well, and the Gil formulation. 
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Figure 13. The Hubble residuals as a function of metallieity converted from mass (left), and the fundamental metallicity relation (right) 
for Gil. The blue line is the linear fit to the Hubble residuals, and the red dots are the median values of the Hubble residuals. The 
likelihood-ratio test shows that the value of the best fit line with the FMR metallicities is not statistically preferred. Although the 
is lower, the significance of the improvement is less than 95% confident. 

for SED fitters and SN la studies alike, especially in light of the potential power of photometric fitting and the FMR, 
an improvement in the estimation of the SFR from photometric SED fitters is of paramount importance. The SFR 
values from SED fitting are highly dependent on the star-formation history used in creating the templates, and the 
constraints on the age of the template that are appropriate for stellar mass estimation are not generally appropriate 
for SFR estimation. Further study of photometric formulation of the FMR (for any purpose) will be highly dependent 
on improvements in the estimation of the SFR. 

4.3. Extending the FMR to Lower Metallicity 

While we have demonstrated the power of photometric fitter s in improving metal licity estimates, and even Hubble 
residuals, our calibration sample was relatively low redshift. .Mannucci et al.l ()2010l ) claims no evolution of the FMR 
out to redshift 2.5, however, the number of galaxies beyond the SDSS sample, which only goes out to reds hift 0.3, was 
small. As redshift increases, the mass decreases and the SFR increases, on average (jHavden et al.|[20T2l ). and so the 
value of the parameter /ifmr decreases as does the metallicity. It would be unwise to apply the FMR to galaxies that do 
not fall within the region of /ifmr and 12-|-Log(0/H) values that were used in the calibration sample. Man nucci et al.l 
1) was able to show that the FMR extends smoothly down to 12+Log(0/H) of aro und 8, usin g the lower-mass 
lower-mctallicity SDSS galaxies that fell into bins that contained less than 50 galaxies. iCresci eFal 
the F MR t o about a redsh i ft of .8, finding no evidence for evolution. This confirms the results in 
(|2010f) and iMannucci et al.l ()2011f) with a higher level of s ignificance 



(f2dn extends 



Mannucci et al.l 



If the FMR is reliable out to redshift 2.5, as implied by IMannucci et all (|2010[ ) and ICresci et al.l (|2012l ). this would 
become a powerful tool for reliably estimating metallicities of a significant fraction of the universe. More work must 
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Figure 14. The likelihood ratio test probability as a function of the FMR parameter a, for Gil. The probability value is a result of using 
the likelihood ratio test statistic A, from equation [5] as the test statistic in the probability distribution with 1 degree of freedom. This 
compares each value of a against the mass alone (shown by a vertical line at a=0). This minimization was performed using the Hubble 
residuals, and is mathematically independent of the FMR. The region of parameter space where the probability is below 5% agrees with 
the best-fit a value of the independently measured FMR (indicated by a star). This indicates the improvement of the FMR in reducing 
the scatter of the Hubble residuals and is evidence that the mass correlation is tracing the metallicity. 

be done to both extend the FMR to larger samples of low metallicity galaxies, and probe potential evolution of the 
FMR extending back in cosmic history. 



5. CONCLUSIONS 

The fundamental metallicity relation shows that a combination of stellar mas s and SFR improves m etallicity estimates 
beyond those derived from stellar masses alone, i.e. the Tremonti relation (Tre monti et al.l 12004"). We have demon- 
strate d the capability of two photometric methods, Z-PE G and minirn i zation of FSPS models from iGupta et al.l 
()2011|) . to reproduce the fundamental metallicity relation of lMannucci et ahl (|2010[ ). and significantly reduce the scatter 
in metallicity estimates. This allows for accurate metallicity estimates based on multi-band photometry alone. 

We have demonstrated the necessity of calibrating the F MR for the photometric fitter that will be used to estimate 
the stellar mass and SFR. In the case of lGupta et al.l (|2011[ ) fits, the relat ively small va r iation in available specific SFR 
values greatly changes the calibration of t he FMR. Using t he Z- PEG or IGupta et al.l (120111 ) galaxy parameters with 
the specific formulation of the FMR from iMannucci et al.l (|2010( ) would give results that were at best systematically 
biased in some way, and at worst nonsensical. 

We have calibrated the FMR for both the Z-PEG and Gil photometric fits, and fit a star-forming subsample of the 
206 SDSS-II SN la host galaxies with these fitters in the SDSS ugriz bands. We have demonstrated that using the 
FMR parameter /.jfmr = log(mass) — ax log(SFR) improves Hubble residual correlation beyond the stellar mass alone 
at a statistical significance of at least 99%, using the likelihood ratio test. We note again that this sample contains only 
relatively high SFR galaxies. Still, the improvement in Hubble residual correlation suggests that the mass correlation 
confirmed by many studies is very likely tracing an underlying metallicity correlation. Generally, metallicity effects 
on SN la peak brightnesses should affect the amount of radioactive Nickel-56 produced, and therefore be corrected 
for by light-curve width corrections. If indeed the mass correlation traces a correlation with the metallicity, then this 
correlation is separate from the light-curve width correction and indicates other physical processes, possibly indicating 
that the Phillips relation itself is metallicity dependent (.Kasen et al. 2009). The calibration of the FMR for photometric 
fitters will allow ever greater precision in the measurement of cosmological distances using SNe la, as well as improving 
the quality of metallicity studies for galaxies where measuring spectra is difficult or nearly impossible. 

Future SN la searches will need to account for the correlation of SN la Hubble residuals with the FMR and metallicity 
illustrated in this work. Where possible, large samples of spectroscopic metallicity values for SN la host galaxies should 
be obtained. Multi-wavelength photometry that allows accurate estimation of the stellar mass and SFR of the host 
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galaxies will be instrumental in using the FMR to correct SN la distances, where spectroscopic study is prohibitively 
expensive, as well as furthering our understanding of the SFR estimates from SED fitting. As SN la studies extend to 
high-redshift, the effect of host galaxy correlations on SN la distances can systematically change, and building on the 
work presented here will be a significant effort for the next generation of SN la surveys. 
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